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I. INTRODUCTION 

In general relativity, non-gravitational processes 
are included via the energy-momentum tensor for 
matter. Various non-gravitational processes de- 
termine, apart from other physical characteristics 
such as radiation, the relation of density and pres- 
sure of matter. The temporal evolution of matter 
is to be determined from such a relation, and from 
other physical relations, if any. 

By physically realizable gravitational collapse, we 
mean gravitational collapse of matter that leads 
matter, step by step, through different "physical" 
stages of evolution, namely, from dust to mat- 
ter with pressure and radiation, to matter with 
exothermic nuclear reactions etc. 

Then, the spacetime of "physically realizable" 
collapse of matter must be able to begin with any 
stage in the chain of evolution of matter under the 
action of its self-gravity. The temporal evolution 
from any "initial" data, any "physical" stage in 
question, is to be obtained from applicable non- 
gravitational properties of matter. This is as per 
the principle of causal development of data. 

In a recent work we obtained a spherically 
symmetric spacetime by considering a metric sep- 
arable in co-moving coordinates and by impos- 
ing a relation of pressure p and density p of the 
barotropic form p = ap where a is a constant. 
Such a relation determined only the temporal met- 
ric functions of that spacetime. 

Therefore, although the temporal metric func- 
tions for this spacetime were determined in |lj us- 
ing the above barotropic equation of state, same 
metric functions arc determinable from any rela- 
tion of pressure and density. (See later.) 

This spacetime admits a hyper-surface orthogo- 
nal, radial, Homothetic Killing Vector (HKV). (See 
later.) Hence, it will be called a radially homoth- 



etic spacetime. Since this spacetime describes ap- 
propriate "physical" stages of evolution of spher- 
ical matter, we have argued in Q that it is the 
spacetime of physically realizable spherical collapse 
of matter. We have also shown in [|| that it is a 
Petrov-type D spacetime. We note that all general 
relativistic black hole spacetimes are Petrov-type 
D spacetimes. 

In H , we studied the shear-free problem for the 
sake of physical understanding of the issues in- 
volved in the study of the spherical gravitational 
collapse in the spacetime of [Q . 

In this paper, we consider a radially homothetic 
spherically symmetric problem for non-singular 
and non-degenerate matter data in its full "physi- 
cal" generality, ie, with shear and radiation. Fur- 
ther, we also obtain the source-free electromag- 
netic fields in this spacetime explicitly. 

We describe, in §y, the spacetime under con- 
sideration and its properties. Next, we summarize 
the Hertz-Debye formalism in 



HI. The electro- 



magnetic fields in _aradially _homothetic spacetime 
are obtained in 



IV 



In §[y|, we discuss the as- 
trophysical relevance and other implications of the 
results obtained in these works. 



II. SPACETIME METRIC 

In general, a HKV captures [Q the notion of the 
scale-invariance of the spacetime. If, in terms of 
the chosen coordinates, a homothetic Killing vec- 
tor X has component only in the direction of one 
coordinate, the Einstein field equations separate 
for that coordinate, generating also an arbitrary 
function of that coordinate. This is the broadest 
(Lie) sense of the scale-invariance leading not only 
to the reduction of the field equations as partial 
differential equations to ordinary differential equa- 
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tions but leading also to their separation. 

A spherically symmetric spacetime has only one 
spatial scale associated with it - the radial distance 
scale. Therefore, for a radially nomothetic space- 
time, the metric admits one arbitrary function of 
the radial coordinate. We then obtain for a ra- 
dially homothetic spacetime arbitrary radial char- 
acteristics for matter. That is, due to the radial 
scale-invariance of the spherical spacetime, matter 
has arbitrary radial properties in a radially homo- 
thetic spacetime. 

In co-moving coordinates, a radially homothetic, 
spherical spacetime admits a spacelike HKV of the 
form 



X a = (0,^,0,0) 



(1) 



and the spacetime metric is given by 

ds 2 = -y 2 dt 2 + 1 2 {y') 2 B 2 dr 2 + y 2 Y 2 dVL 2 (2) 

with a prime indicating a derivative with respect 
to r, B = B(t), Y = Y(t) and 7 being a constant. 
(We absorb the temporal function in g u by redefi- 
nition of the time coordinate.) 
The Ricci scalar for (0) is: 



K = 



AYB IB 



G 



y 2 YB y 2 B y 2 1 2 B 2 

2 2Y 2 AY 



yiy 2 y 2 Y 2 y 2 Y 



(3) 



The non-vanishing components of the Weyl ten- 
sor for (0) are: 



C, 



a 



LOIO 



a 



t<jit<p 



GrdrB 



3 

y 2 Y 2 

6 



F(t) 



F(t) 



sin aUtete 
B 2 j 2 Y 2 (y') 2 
(3 



F(t) 



C r( p r( p — sin 2 9 C r g r g 
y 2 Y 4 sin 2 1 

0< 



(4) 

(5) 
(6) 

(7) 
(8) 



F(t) (9) 



where 



, N Y Y 2 1 B BY 
F ^ = Y-^-Y^-B + BY (10) 

In what follows, we shall assume, unless stated 
explicitly, that there are no singular initial- data 
and that there are no degenerate situations for the 
metric See later for the singularities and the 
degeneracies of the metric ((2]). 



It is well-known that spherical spacetimes are 
either Petrov-type O or Petrov-type D. Type-D 
spacetimes are not conformally flat. As can be 
easily verified, the spacetime of (||) is not con- 
formally flat for non-singular and non-degenerate 
data. That the spacetime of (^) is of Petrov- 
type D can also be seen differently by verifying || 
that only $2 = — C a b c d^ a, m h rh c n d is non- vanishing 
for the metric ^ where I, n, m and m are the 
Newman- Penrose || tetrad vectors for it. 



A. Elementary flatness and center 

Now, the spacetime of (||) is required, by defini- 
tion, to be locally flat at all of its points including 
the center. 

In the case of (||), a small circle of coordinate 
radius e with center at the origin has circumference 
of 2ire. On the other hand, the circle has the proper 
radius 7 y' e. Then, requiring that the ratio of the 
circumference to the proper radius of the circle to 
be 2ir in the neighborhood of the origin, we obtain 
the condition for the center to possess a locally flat 
neighborhood as 



y 



l/ 7 



(11) 



This condition must be imposed on any y(r). With 
this condition, (|ll|), the HKV of metric (||) is, at 
the center, y\ r =od/d r . 

Now, y(r) is the "area radius" in (||). When 
y\ r =o 7^ 0, the orbits of the rotation group 50(3) 
do not shrink to zero radius at the center for (g). 
Consequently, the center is not regular for (0) when 
y\ r =o 7^ although the curvature invariants remain 
finite at the center. 

Also, when y\ r =o = 0, the center is regular for 
the spacetime of (|2]). But, the curvature invariants 
blow up at the center, then. 

It is well-known |^| that the center and the initial 
data for matter, both, are not simultaneously reg- 
ular for a spherical spacetime with hyper-surface 
orthogonal HKV. Therefore, the spacetime of (||) 
does not possess a regular center and regular mat- 
ter data, simultaneously. 

However, the lack of regularity of the center of 
(E) for non-singular matter data is understandable 
since the orbits of the rotation group do not 
shrink to zero radius for every observer. It is a 
relative conception and the co-moving observer of 
(||) is not expected to observe the orbits shrink to 
zero radius. (See also later.) 
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B. Singularities of spacetime 



via a non-singular coordinate transformation 



Clearly, we may use the function y(r) in (Q) as 
a new radial coordinate - the area coordinate - as 
long as y' ^ 0. However, the situation of y' = 
represents a coordinate singularity that is similar 
to, for example, the one on the surface of a unit 
sphere where the analogue of y is sin6> Q|. The 
curvature invariants do not blow up at locations 
for which y' = 0. 

The genuine spacetime singularities of the strong 
curvature, shell-focussing type exist when either 
y(r) = for some r or when the temporal func- 
tions vanish for some t — t s . 

There are, therefore, two types of curvature sin- 
gularities of the spacetime of @, namely, the 
first type for B(t s ) — and, the second type for 
y{r) = for some r. 

Note that the "physical" radial distance corre- 
sponding to the "coordinate" radial distance 6r is 



I = 1 {y')B5r 



(12) 



Then, collapsing matter forms the spacetime sin- 
gularity in (Q) when B(t) = is reached for it at 
some t — t s . Therefore, the singularity of first 
type is a singular hyper-surface for (^). 

The singularity of the second type is a singular 
sphere of coordinate radius r. The singular sphere 
reduces to a singular point for r = that is the 
center of symmetry. For y(r) =0 for some range 
of r, there is a singular thick shell. Singularities 
of the second type constitute a part of the initial 
data, singular data, for the evolution. 



C. Degeneracies of the metric (g) 

The metric (||) has evident degeneracies when 
y(r) = 0, y{r) — > oo either on a degenerate sphere 
of coordinate radius r, for some "thick shell" or 
globally. The degeneracy y(r) = is equivalent to 
an infinite density while the degeneracy y(r) = oo 
is equivalent to vacuum. (Sec later for the expres- 
sion of density of matter in the spacetime of (g).) 
Another degeneracy occurs for y(r) = constant for 
some "thick" shell or globally. This degeneracy 
corresponds to uniform density. 



T 



l(t) exp 



k(t) exp 



F^dr 



(14) 



F^dr 



We note that, if we invoke (Q) for (§), the re- 
sulting metric will not be diagonal. The imposition 
of diagonality of the metric will require a relation- 
ship between l(t) and k(t). Such a relation can, of 
course, always be imposed. 

Hence, the metric (||) can always be trans- 
formed, under non-singular coordinate transforma- 
tions (0), to a form which admits a HKV ( |l3"| ) 
in the transformed coordinates. The transformed 
metric under consideration is, therefore, 

ds 2 = -P 2 dT 2 + Q 2 dS 2 + S 2 Z 2 dil 2 (15) 

where P,Q,Z are the metric functions of the self- 
similarity variable T/S or S/T. Note that for the 
transformed metric the radiation or the heat flux 
is, in general, non-vanishing. 

For (|l5|), we are therefore led to consider the 
spacetime singularity at S — and T = 0. 

As we noted earlier, a relationship exists be- 
tween l(t) and k(t) of (Q). Thus, S = and 
T = for @ corresponds to l(t) = for (§ 
when l(t) oc k(t) and l{t) = 0. There is thus no 
constraint on the radial function yir) in (|^) that 
it should vanish. Consequently, y(r) ^ at r = 
or, for that matter, at any r, is permissible. Then, 
y\ r =o is arbitrary. 

However, 5 = and T = also corresponds 
to y(r) = in ( |i"4| ) with l(t) not proportional to 
k(t). However, assuming y(r) = leads to ini- 
tially singular density. Such a curvature singu- 
larity will then always exist on any hyper-surface 
t = constant in the spacetime of (||). 

We note that we may begin with the metric form 
( p^l ) even when the transformations (Q) are singu- 
lar. The metric (|l^) is then not reducible to the 
form (0). For such a spacetime, the Einstein field 
equations reduce to ordinary differential equations. 
However, the resultant equations are not entirely 
separated in terms of the variables T and S. 



D. Self-similarity Nature 

Any vector (|l|) can always be transformed 
into 



X a = (T, 5,0,0) 



(13) 



E. Mass function 

The mass function for the spacetime of (|2|) can 
be defined as 



(16) 
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where m(r, t) denotes the effective total energy per 
unit mass of a fluid element labelled by the co- 
moving radial coordinate r at co-moving time t. in 
the spacetime. Note that it includes the "effective" 
contribution due to the flux of radiation or heat in 
the spacetime. 

Note that, using (|l^), we may rewrite the metric 
(0) as: 



ds 2 = - dr 2 



R' 2 dr 2 



1 - {2m/ R) 



R 2 dn 2 



(17) 



where R = R(r,r) — yY, m = m(r, r) and 
dr = ydt. Now, the metric (g) is recognizable as a 
generalization of the Tolman-Bondi dust metric to 
include pressure and radiation. 

Now, we note the following. An asymptotic ob- 
server of the Schwarzschild spacetime does not see 
a sphere, drawn around the central mass point, 
shrink to a zero radius as a result of the red-shift 
of the sphere becoming infinite at the gravitational 
radius. Consequently, the "area radius" of the 
sphere does not shrink to a zero for the asymp- 
totic observer but it appears to this observer that 
the minimum area radius of the sphere is r = 2M, 
where M is the central Schwarzschild mass. For 
the asymptotic observer, r = 2M is then the "cen- 
ter" of the spacetime. 

A co-moving observer, as a cosmological ob- 
server in (Q), is "equivalent" to the asymptotic 
observer of the Schwarzschild geometry. Hence, 
for the co-moving observer of (|2|), the area of any 
sphere cannot shrink to zero radius at the center 
since there is always an equivalent mass point at 
the center for that observer. (See also Q.) 

Therefore, the lack of regularity of the center of 

tfor non-singular matter data is understandable 
since the orbits of the rotation group do not 
shrink to zero radius for the co-moving observer of 
(2). Note also that the singular hyper-surface of 
(2) is then the infinite-mass singularity. 



F. Field equations 

For completeness, we reproduce the arguments 
of § here. 

Now, define the co-moving time-derivative 



D t = U a 



d Id 
dx a y dt 



where 



y 



(18) 



(19) 



is the four-velocity of the co-moving observer. 
Then, the radial velocity of fluid with respect to 



the co-moving observer is 

V r = D t (yY) = Y 



(20) 



where an overhead dot has been used to denote a 
time derivative. 

The co-moving observer is accelerating for (||) 
since 



u a = u a . b u b 



o, -, o, 

V 



(21) 



is, in general, non-vanishing for y' ^ 0. The ex- 
pansion is 



(22) 




The Einstein tensor for (|2|) is: 



G tt 



1 1 

Y 2 ~ <y 2 B 2 



Y 2 



'BY 



(23) 



y- 



-2- - — 

Y ~ Y 



1 

Y 2 



Gee 

Gifted 

Gtr 



,B 



= -YY — Y — - Y ■ 



B 



YB 
~B~ 



Y 2 



j 2 B 2 



sin 2 G e 
By 



(24) 

(25) 
(26) 
(27) 



Now, to see, and only to see, that B is related 
to the flux of radiation in the spacetime of (|J), 
we may consider || matter to be described by the 
energy-momentum tensor 



fab = (p + p)U a U b + pg ab 



(28) 



and consider that it fails to satisfy a local con- 
servation law due to the emission of radiation that 
escapes radially along the radial null vector l a . For 
the radiation, we may then assume the "geometri- 
cal optics" form 



E ab = Qgap 



(29) 



with Q being the energy density of radiation or the 
energy flux density in the rest frame of the fluid. 
Then, it is seen that Q oc B. Thus, the radiation- 
flux depends on B. 

Now, define the quantity 



a= „\=a 2 2 = --a\ = — (---] (30) 
2 3y \ Y B / V ' 
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Here, cr a b represents the shear-tensor of the fluid 
and the shear-scalar is given by V6 a. The space- 
time of (Q) is shearing when B(t) is not propor- 
tional to Y(t). Therefore, the spacetime of (||) is, 
in general, shearing and radiating, both. 

Now, we note that Penrose is led to the Weyl 
hypothesis on the basis of thermodynamical con- 
siderations, in particular, those related to the ther- 
modynamic arrow of time. On the basis of these 
considerations, we may consider the Weyl tensor 
to be "some" sort of measure of the entropy in the 
spacetime at any given epoch. 

Then, for non-singular and non-degenerate data 
in (||), the Weyl tensor of (||) blows up at the sin- 
gular hyper-surface of (||) but is "constant" at the 
"initial" hyper-surface since Y = B = for the 
"initial" hyper-surface |2j. 

This behavior of the Weyl tensor of (||) is in con- 
formity with Penrose's Weyl curvature hypothesis 
@. Thus, the spacetime of @ has the "right" 
kind of thermodynamic arrow of time in it. 



Stages of collapse 

Now, we turn to steps of collapse of mat- 
ter in the spacetime of (||) for non-singular and 
non-degenerate data. For non-singular and non- 
degenerate data, we then have a "cosmological" 
situation - continued spherical collapse of matter 
from the assumed "initial" state. 



Step I - Evolution of dust 

We begin with collision-less and pressureless 
dust matter with density distribution given such 
that y(r) > 1 everywhere on the initial hyper- 
surface in the spacetime of (||). Emission of ra- 
diation and, hence, radiation itself is not expected 
in such dust. 

Then, self-gravity leads to mass or energy-flux 
in the radial direction. But, this is not the flux 
of radiation. Therefore, there is no mass-flux in 
the rest frame of collapsing dusty matter, but it is 
present for other observers in the spacetime. 

Then, for vanishing flux of radiation, we have 
from (|2j), B = constant = B Q . 

Then, the co-moving density, p, of dust is 



1 

2 

y 



Y 2 



1 

y2 



1 



1 2 B 2 



(31) 



and the function Y(t) is determined by the condi- 
tion of vanishing of the isotropic pressure: 



Here, £ = 5/"f 2 B 2 , a positive constant. 

A solution of this equation is obtainable as 



dY 



s/-l + C/5Y 2 + CoY- 1 / 2 



t-t 



(33) 



where c D is constant. Since Y is the radial velocity 
of matter for the co-moving observer, we require 
that solution for which Y — >• for t — > — oo. 

That the dust exists in the spacetime of (0) is 
not surprising since it is a generalization of the 
Tolman-Bondi dust metric. 



Step II - Evolution with pressure and radiation 

Next stage of collapse is reached when particles 
of dust begin to collide with each other. Negligi- 
ble amount of radiation, but existing nonetheless, 
is expected from whatever atomic excitations or 
from whatever free electrons get created in atomic 
collisions in such matter. Therefore, dusty matter 
evolves into matter with pressure and radiation, 
both simultaneously non- vanishing. 

The energy-flux can no longer be removed by 
going to the rest frame of matter. 

Now, pressure and radiation, both, get simulta- 
neously switched on in the spacetime of (||) when 
B(t) 7^ 0. This is as per the expectation that 
dusty matter evolves to one with simultaneous oc- 
currence of pressure and radiation, both. 

Then, the co-moving density, p and isotropic 
pressure, p, are given by 



1 



Y 2 



3 7 2 B 2 



YB 
YB * 

1 

~ W 2 

2B 
W ~ 



1 

Y 2 ~ 

AY 
~ SY 

Y 2 

3Y 2 ' 



1 



2YB 
3YB 



(34) 



(35) 



From (|3) and (p|), we obtain 



Y 
Y 



B 
B 



i 2 B 2 



-(p + 3p) (36) 



AYY + Y 2 + 1 - C y 







(32) 



Then, from (|3q), the relation of pressure and 
density of matter is the required additional "phys- 
ical" information. Also required is other relevant 
"physical" information to determine the radiation 
generation in the spacetime of (||). This is a non- 
trivial task in general relativity just as it is for 
Newtonian gravity. 

Note that the radiation may be "negligible" but 
what is important to considerations here is its pres- 
ence in the spacetime. 



G 



Step III - Stellar object 

At some further stage of evolution, radiation 
from "central part" of collapsing matter may 
become non-negligible and the temperatures in 
the central region may become appropriate for 
exothermic, thermonuclear reactions. With the 
onset of exothermic thermonuclear reactions, a 
"shinning" star is born in the spacetime. 

The exothermic thermonuclear reactions in the 
stellar core may support the overlying "stellar lay- 
ers" and such a stellar object may "appear" grav- 
itationally stable. 

But, the spacetime continues to be dynamic 
since radiation is present in it. The central stellar 
object may also accrete matter from its surround- 
ing while emitting radiation. 

Now, as and when "heating" of the overlying 
stellar layers decreases due to changes in exother- 
mic thermonuclear processes in the core of the star, 
the self-gravity of the stellar object leads to its 
gravitational contraction. These are, in general, 
very slow and involved processes. 

Gravitational contraction leads to generation of 
pressure by compression and by the occurrence 
of exothermic thermonuclear reactions involving 
heavier nuclei. The star may stabilize once more. 

This chain, of gravitational contraction of star, 
followed by pressure increase, followed by subse- 
quent stellar stabilization, continues as long as 
thermonuclear processes produce enough heat to 
support the overlying stellar layers. 

The theory of the atomic nucleus shows that 
exothermic nuclear processes do not occur when 
Iron nucleus forms. With time, the rate of heat 
generation in iron-dominated-core becomes insuf- 
ficient to support the overlying stellar layers which 
may then bounce off the iron-core resulting into a 
stellar explosion, a supernova. 

Then, many, different such, stages of evolution 
are the results of physical processes that are un- 
related to the phenomenon of gravitation. These 
are, for example, collisions of particles of matter, 
electromagnetic and other forces between atomic 
or sub-atomic constituents of matter etc. 

As an example, let some non-gravitational pro- 
cess, opposing collapse, result into pressure that 
does not appreciably rise in response to small con- 
traction of the stellar matter. That is, pressure 
does not appreciably rise when gravitational field is 
increased by a small amount Then, the collapse of 
a sufficiently massive object would not be halted by 
that particular non-gravitational process. There- 
fore, a mass limit is obtained in this situation. For 
example, electron degeneracy pressure leads to the 
Chandrasekhar limit [^TJ. 



Clearly, some of the non-gravitational processes 
determine the gravitational stability of physical 
objects. This is true in Newtonian gravity as well 
as in general relativity, both. 

If we consider that density alone does not de- 
termine pressure completely but that the isotropic 
pressure is a function of density and temperature, 
both, then we need another equation, from thermo- 
dynamical considerations, perhaps. In this case, 
the relation of pressure and density can be consid- 
ered to be a function of time. 

But, since the relation of pressure and density is 
arbitrary for (|J) , a changing pressure-density rela- 
tionship is allowed for it, we note. To provide for 
the required physical information is, once again, a 
non-trivial task. 



Evidently, to provide for the required informa- 
tion of "physical" nature is a non-trivial task in 
general relativity just as it is for Newtonian grav- 
ity Q. The details of these considerations are, of 
course, very involved and have been left out of the 
considerations of the present paper. 

We note that, in general, the relation of pres- 
sure and density at extremely high densities is not 
known. However, it can be surmised that, in the fi- 
nal stages of collapse, the collapsing matter will be 
ultra-relativistic and will end up in the singularity 
as such. Then, the relation of pressure and density 
of such matter may be expected to remain fairly 
unchanged in the final stages of the gravitational 
collapse. That is to say, the spacetime of collaps- 
ing matter should describe a relation of pressure 
and density of ultra-relativistic matter that is not 
changing in the final collapse stages. 

Consequently, it seems reasonable to treat the 
thermodynamic state of collapsing matter near 
the spacetime singularity by a relation of the 
"barotropic" form p = ap where a is a constant 
characteristic of the collapsing matter. Therefore, 
by assuming this equation of state, we may obtain 
the temporal metric functions in (|J) to study the 

ftvitational collapse in its final stages. But, for 
, such a possibility of final stage of collapse is 
approached only asym ptotica lly, that is for infinite 
co-moving time. (See § IIH 1 later.) 

But, it is clear that the field equations deter- 
mine only the temporal functions of (|[) from any 
suitable energy-momentum tensor including that 
of electromagnetic fields, if any. 

The temporal functions B(t) and Y(t) are to be 
obtained from the properties of matter such as a 
relation of pressure and density, the rate of loss of 
internal energy to radiation, processes of quantum 
mechanical nature etc. 

Moreover, it is also clear that matter will con- 
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tinue to pile up on such a star in a "cosmolog- 
ical setting" and, hence, such a star will always 
be taken over any mass-limit in operation at any 
stage of its evolution. The evolution of collapsing 
matter will, asymptotically, lead to the singular 
hyper-surface of the spacetime of (|^) . 

The radial dependence of matter properties is 
"specified" as 1/y 2 but the field equations of gen- 
eral relativity do not determine the metric function 
y(r) in (g). 

Therefore, the radial distribution of matter is ar- 
bitrary in terms of the co-moving radial coordinate 
r. This is the "maximal" physical freedom com- 
patible with the assumption of spherical symme- 
try, we may note. Note, however, that the phys- 
ical generality here is not be taken to mean the 
"geometrical" generality. 



G. Absence of null or one-way membrane 

A spherical surface r = constant in the geometry 
of (|^) has a normal vector 



n a = (0,l,0,0) n a n a 



1 



(37) 



Within the range (0, oo) of the co-moving radial 
coordinate r, the character of n a does not change 
from spacelike to null to timelike in the spacetime 
of (||). Then, the norm n a n a does not vanish at 
any r in non-singular and non-degenerate cases. 

To see the same differently, the coordinate speed 
of light in the spacetime of (0) is 



dr 



= ± 



7 (y')B 



(38) 



This speed cannot vanish for non-singular and non- 
degenerate cases. At the singular hyper-surface, 
the coordinate speed of light becomes infinite. 

Therefore, in the spacetime of (|J), there can- 
not be a spherical, spatially finite, null membrane 
or a one-way membrane, ie, a black hole in the 
usual sense of the term, in non-singular and non- 
degenerate cases. 



H. Spherical collapse, singular hyper-surface 
and the infinite red-shift surface 

We emphasize that we are using non-singular 
and non-degenerate data for (||) . A co- moving ob- 
server in (H) is then "a cosmological observer" . The 
four-velocity of matter fluid with respect to a co- 
moving observer is: 



0,0) 



v r = 



(39) 



We then obtain from the metric (||): 
1 



(1, V"\0,0) 



(40) 



A = 1 - 7 2 (^j B 2 {V r f (41) 

Now, if dr CM is a small time duration for a co- 
moving observer and if dr RF is the corresponding 
time duration for the observer in the rest frame of 
matter, then we have 



dr r 



dr Rf 



From (f42|), we also get the red-shift formula 

V C m = Vrf %/A 1 + 2 = 1 



Vcm 



(42) 



(43) 



in the spacetime of (Q) where v CM is frequency of 
a photon in the co-moving frame, v RF is the fre- 
quency in the rest frame and z is the red-shift of 
the photon. Then, A = is the infinite red-shift 
surface that is, however, not a null membrane. A 
co-moving observer waits for an infinite period of 
its time to receive a signal from the rest-frame ob- 
server when A = 0. 

Then, we distinguish regions of the spacetime of 
©as 

(A > 0) \ 7 (y')BV r \ < y (44) 
(A = 0) h(y')BV r \ = y (45) 
(A < 0) h(y')BV r \ > y (46) 

Now, the geodesic equations of motion for (||) 
are easily obtainable. The r-equation is: 



Ts (^yy'r] 



2C 



(47) 



where an overhead tilde denotes derivative with 
respect to the affine parameter s. Then, 



2Cs + ki 
7 2 B 2 y y' 



(48) 



where k\ is a constant of integration. 

For the motion of the particle in the equatorial 
plane 9 = 7r/2, the solution of the t-equation of 
motion, using (J4q) in the lagrangian of motion, is 



t 




1 



{2Cs + k 1 ) 



Then, since V r = dr/ydt, we have 
(fci - s) 2 



(49) 



V r V r = 



y 2 j 2 B 2 + (fci - s) 2 



(2£ = -1) (50) 



where s is the affine parameter along the geodesic 
and we have V r = at s = k\. 



s 



1. Singular hyper-surface 

Clearly, V r V r = 1 for B(t B ) = at t = t s , ie, the 
velocity of the particle with respect to a co-moving 
observer is the speed of light at the singular hyper- 
surface of (§) - B(t s ) = at t = t s . This is, of 
course, happening only asymptotically. 

But, from @, A = 1 for B(t s ) = 0. A co- 
moving observer also moves with the speed of light 
at the singular hyper-surface. After all, matter ev- 
erywhere should become relativistic as the central 
mass condensate continues to grow (due to accre- 
tion onto it) to influence the entire spacetime to be- 
come relativistic everywhere. This too is happen- 
ing only asymptotically. But, the singular hyper- 
surface is not the infinite red-shift surface. 



2. Infinite red-shift surface 

A spatially bounded infinite red-shift surface 
"occurs" in (H) when A = 0. Since 



7 yg 2 + (fci- g ) 2 (i-^) 

j 2 y 2 B 2 + (ki - s) 2 



(51) 



this requires y(r) < 1. 

Thus, the initial "density distribution" , from 
y(r), decides whether (||) has an infinite red-shift 
surface or not, ie, whether y 2 — V r V r . This is 
possible for V r V r « 1, but for sufficient mass 
concentration at the center, ie, for y < 1. 

But, a co- moving observer does not see the in- 
finite red-shift surface form. Thus, an initially 
"small" matter density (y > 1) cannot, for a co- 
moving observer, become "large" enough that the 
infinite red-shift surface forms (y < 1). 

Matter piles up at the center of the spacetime. 
But, the density at the center does not become 
infinite at any finite co-moving time. 

Now, let y < L ie, let there be an infinite red- 
shift surface in (||) at the initial epoch itself. (For 
y\ r —o = 0, there is an infinite-density singularity 
at the center.) Then, from (48), it follows that a 
particle can reach and cross this infinite red-shift 
surface along its radial geodesic. 

But, such a central region - "y < 1 apparition" 
- cannot communicate to its exterior since the co- 
moving time duration is imaginary in this region. 
Therefore, the interior of an infinite red-shift sur- 
face is causally disconnected from the rest of the 
spacetime at the initial time itself. 

Then, matter external to y < 1 region eventually 
collapses onto this central "apparition" without a 
co-moving observer seeing matter enter it. The 
spacetime of (H) then describes the accretion of 



matter onto the central "apparition" - the spatially 
bounded infinite red-shift surface. 

We shall refer to the infinite red-shift surface, as 
described above, as the black hole surface and its 
interior as an interior of a black hole. This is the 
conception of a black hole that arises in the space- 
time of (^) for non-singular and non-degenerate 
data. But, here, a black hole is not a null mem- 
brane or horizon. 



Apparent horizon 



A radially outgoing null vector of (0) is 



e°d =-- 

ydt 



1 d 

^yy'B dr 



(52) 



Light gets trapped inside a particular radial coor- 
dinate r when the expansion of the above principle 
null vector vanishes at r. The formation of the out- 
ermost light-trapping surface or the apparent hori- 
zon is then obtained by setting the expansion of 
( [52"1) to zero. 

The zero-expansion of ( |52| ) yields a condition 
only on the temporal metric functions as 



B Y _ _3_ 

B + Y ~^B 



(53) 



The condition ( j53| ) implies an "instant of time" . 

This is seen as as follows. An outgoing photon 
moves along the trajectory 



dr 
~dt 



y 



i{y')B 



and crosses a sphere of coordinate radius r at co- 
moving time t. 

The mass inside this sphere is given by (|l6|). 
Now, the equation of the light trajectory can be 
used to express the mass function m as a function 
of either r alone or t alone. The light trapping 
mass is then obtained by setting 2m = yY. Then, 
for every value of r there is some t and vice versa 
for which 2m = yY . Thus, condition ( |53| ) implies, 
in essence, an "instant of light trapping" . Essen- 
tially, for non-singular and non-degenerate data in 
(^|), every instant, of the co- moving time, is an in- 
stant of light trapping. 

Alternatively, let a spherical light front be emit- 
ted from the center of symmetry As it travels radi- 
ally outwards, it brings in more mass to its interior. 
When sufficient mass is in the interior, light trap- 
ping occurs. For non-singular and non-degenerate 
data in (Q), we can always draw a sphere contain- 
ing enough mass that can trap light. This is the 
essence of the statement that "Every co-moving 
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instant is a 'Light Trapping Instant' in this space- 
time" . In a sense, every observer is inside some 
light trapped sphere in (|J) for non-singular and 
non-degenerate data. 



I. Shell black hole 

Now, an interesting possibility is that of (||) con- 
taining many concentric infinite red-shift surfaces. 
Then, what we have here is the possibility of shells 
of black holes! 

Intuitively, this is the only possibility that can 
arise in spherical symmetry apart from that of a 
single spherical black hole as considered earlier. It 
is also clear that this possibility arises only as "ini- 
tial data" in the spacetime of (|^). 

To analyze such a "shell" black hole, we will, of 
course, be required to use the gaussian coordinate 
system since y' = at some radial locations. We 
recall that y' = is a coordinate singularity of the 
metric (|J). 

That the shell black hole is obtained in the 
spacetime of (^) is not a coincidence since (||) is 
the spherically symmetric spacetime with maximal 
"physical" freedom. 



J. When is the singularity of (JiJ) a locally or 
globally naked singularity? 



The radial null geodesic in (g) satisfies 

dt y' 
— = ±7— B 
ar y 



The above differential equation of the radial null 
geodesic does not possess a singularity for non- 
singular and non-degenerate data. 

For non-singular and non-degenerate data, the 
spacetime singularity develops as a result of only 
the temporal evolution of matter in the spacetime 
to the future, ie, as B(t) — > 0. 

But, then the above tangent is vanishing at the 
singularity indicating that there does not exist a 
future-directed, timelike tangent at the spacetime 
singularity. Hence, the curvature singularity of (^) 
is neither globally nor locally visible for the non- 
singular and non- degenerate data. It is a singular 
hyper-surface occurring to the future of every ob- 
server in the spacetime of (^|). 

On the other hand, for singular data, ie, when 
y(r) = for some r, the equation of radial null 
geodesic ( |38| ) has a singularity at that r. 

Then, whether such a singularity or singular 
sphere is visible to any observer or not depends 
on the limit of the quantity y' jy as we approach 



the singular point for B(t) ^ 0. For some func- 
tions y(r) such a limit can be positive making the 
singularity in question a locally naked or a locally 
visible one. 

Therefore, for the spacetime of the visibility 
of the spacetime singularity is determined princi- 
pally by whether we assume the existence of a visi- 
ble singularity at the initial time or not. The tem- 
poral evolution of non-singular and non-degenerate 
data does not lead to a visible singularity in this 
spacetime, we note. 

A short comment on the the possible naked sin- 
gularity solutions of general relativity will not be 
out of place here. 

For non-singular and non-degenerate data in (|J) , 
no naked singularities arise. A black hole, an infi- 
nite red-shift surface, can exist in (|^) as a part of 
this initial data. 

Now, naked singularities can "arise" for singular 
and degenerate data in (|2|). We emphasize that 
for singular and degenerate data we have to aban- 
don the metric (|J) and seek an entirely different 
solution of the field equations of general relativity. 
In cases wherein the HKV of these spacetimes of 
naked singularities is ( filf ) and it is not reducible 
to (Q), the rationale of "pure" radial self-similarity 
is lost for these spacetimes. 

In the absence of "pure" radial scale-invariance, 
it is then not clear what is the "physical" signifi- 
cance of the existence of the HKV (|l3|). Perhaps, 
there is none. 



III. HERTZ-DEBYE FORMALISM 

Having presented the properties of the spacetime 
of (||), we now turn to the Hertz-Debye formalism 
before embarking upon the nature of the electro- 
magnetic fields in the spacetime of (||). 



A. Hertz-Debye potentials 



In flat space, Hertz |12| introduced two vector 
potentials Pg and Pm related to the standard elec- 
tromagnetic potentials $ and A as 



$ = - V • P, 



E 



A = 



dP E 
dt 



V x P M (54) 



and, hence, the bi-vector or the anti-symmetric 
second-rank tensor potential is related by second 
derivatives to the physical fields. 

The gauge freedom associated with the Hertz 
potentials is such as to preserve the source-free 
character of the Maxwell equations while the gauge 
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terms appear as sources in the equations: 

dC 

fe = VxG Qm = — — — Vg (55) 

at 

and 

dW - - - 

R E = - — - Vw R M = -VxW (56) 
at 

where (G, g) and (W, w) are arbitrary 4-vectors. 

This gauge freedom can be used |ll| to reduce 
the Hertz bi-vector to purely radial vectors of the 
form Pe = Pe f and Pm = Pm ? where f is the 
unit radial vector. The functions Pe and Pm 
are the Debye potentials [|l4| and obey, both, a 
wave equation. It should be noted that only the 
"monopole" field is missing in this scheme . 

In essence, the arbitrary, source-free electromag- 
netic field is specified by two scalar functions which 
obey a single, separable second-order wave equa- 
tion. Therefore, a remarkable economy is achieved 
by the Debye potentials. In |L5|, this is expressed 
as: "since a zero rest-mass field possesses two de- 
grees of freedom, no more economical representa- 
tion of the Maxwell field is possible" than that pro- 
vided by the Debye potentials. 

B. Cohen & Kegeles generalization 

Differential forms generalize the flat space 
Maxwell equations to any curved spacetime in a 
natural way. Define the Maxwell 2-form as 

f=\f ab u a Au b (57) 

where / & is the Maxwell tensor and ui a are the 
basis forms. The Maxwell equations are simply 

df = 5f = (58) 

where d is the exterior derivative and 5 = *d* 
is the co-derivative. Here, * is the Hodge dual 
operation. 

The Hertz bi-vector P (2-form) is related to the 
electro-magnetic four-potential 1-form A and the 
Maxwell 2-form as 

A = SP f = dSP = - 5dP (59) 

Then, the equality of the last two expressions in 
( p9| ) requires that 

AP = (dS + 5d)P = (d*d* + *d*d)P = (60) 
where A is the harmonic operator. 



The 2-form gauge terms are: 

Q = dG R = *dW (61) 

where G and W are arbitrary 1-forms. Therefore, 
the wave equation with the gauge terms is 

AP = dG + *dW (62) 

so that the transformed fields are 

/ = d5P -dG^-kdW -SdP (63) 

The transformed fields still obey the source-free 
Maxwell equations as a consequence of the impor- 
tant identities: 

d 2 = S 2 = (64) 

resulting to df = and 5f = 0. 

Equations ([32|) and ( |63| ) provide an elegant and 
fully covariant generalization of the Hertz potential 
formalism to curved spacetimes. 

The problem now consists of determining spe- 
cial bi-vector directions in the spacetime so that 
( |62| ) yields decoupled wave equations for the cor- 
responding components of the potential for some 
choice of the gauge terms (^Tj) . 

In a class of spacetimes, the principal directions 
of the Weyl tensor ||[ [l(| provide such special 
bi- vectors Such special bi-vector directions 

are defined geometrically and independently of the 
Maxwell fields to be computed. In essence, one 
chooses a null tetrad (the Carter tetrad) with one 
null vector aligned along the repeated principal 
null direction of the Weyl tensor of an algebraically 
special spacetime in this scheme. 

This completes our overview of the Hertz-Dcbyc 
formalism or the Cohen-Kegeles formalism p^[ . 
We now turn to the problem of electromagnetic 
fields in the spacetime of (g). 

IV. ELECTROMAGNETIC FIELDS IN 
THE SPACETIME OF ^ 

Choosing an orthonormal tetrad as a basis, the 
spacetime metric (^]) is 

ds 2 = - (uj ) 2 + (u 1 ) 2 + (uj 2 ) 2 + (to 3 ) 2 (65) 

where 

oj° = ydt uj 1 = yy'Bdr (66) 

uo 2 = yYdQ lj 3 = yY sin 8dcj) (67) 

The Hertzian potential 2-form is chosen to be 

P = P B w l) Aw 1 +P M w 2 Aw 3 (68) 
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Under the Hertz-Debye formalism, we obtain a 
wave equation for Pg and Pm, each. 

The electric components are obtained by setting 
Pm = in (|6|) and solving for Pe while the mag- 
netic components are obtained by setting Pe = 
in ( |68| ) and solving for Pm- Since the resultant 
wave equation is identical in both these cases, we 
shall adopt the generic notation \P for the Debye 
potentials Pe and Pm, both |l5| . 

The physical correspondence with the fields is: 



E, 
B 2 



fiO 

hi 



Pi — hz 

B3 = fl2 



(69) 



where the index i ranges from 1,2,3, Ei are the 
electric field components and Bi are the magnetic 
field components. 



Choose (Sec Appendix - Al, for computational 
details) 

* = T en (t)R n (r)P e (cose)e m ^ (70) 



where P^(cos#) is an associated Legendrc function, 
the temporal function T^ n (t) satisfies (A8) and the 
radial function R n (r) satisfies (|A9|). As usual, £ 
and 771 are integers: m to ensure single- valued na- 
ture of e lm ^ and t to ensure that the associated 
Legendre functions do not diverge for cos# = ±1. 
That is to say, the associated Legendre functions 
are polynomials. 



When y' ^ 0, the radial equation ( |A9| ) can be 
written as: 



1 dR n 



d R n 



dy 2 y dy 



:Rn = 



(71) 



and it is an Euler equation. The solutions of this 
Euler equation, for y > 0, are: 



n > R n = ci cos (-\/«lny) + C2 sin (^Jnlny) 

n<0 R n = c 3 y^ + c 4 y~ ^ 

n = R n = c 5 In y + c 6 



(72) 
(73) 
(74) 



where c±, C2, C3, C4, C5 and cq are constants. 

The temporal equation (A£) is of the Fuchsian form and is amenable to series solutions as per the 
theorem of Fuchs. Alternatively, a substitution 



T en (t) 



T tn (t) 



D 



can be used to recast (A8) into the form 



r - ( 1B2 
ln ' Up 2 y 2 



7 2 P 2 2 P 



- H I T en = -W(t)T en 



(n c - n) 



j2 B 2 



(75) 



(76) 



At any given co-moving time, the temporal function W(t) has an inflexion point W(t) = at a critical 
value denoted by n c (£,t) with 



7 2 P 2 
y 2 



n c (£,t) = 



Y 2 B 2 
4P 2 



f 2 p 

2P 



l(£+l) 



(77) 



r 



For n c > n, the function Ti n {t) displays an oscil- 
latory behavior. On the other hand, for n c < n, 
the function Ti n (t) displays an exponential behav- 
ior. At the critical value, n = n r , the function 



Tin oc ((f + i)/VB where C and t are constants. 

In terms of the solutions Te n {t), R n (r), Pe{cos( 
and e 4Tn( ^, the electric multi-pole fields, from ( p3 
are: 



J 



= e -^rT en (t)R n (r)p e (cosey m * 

y z Y z 



Pi = 



12 



| rj y jyt pCf ft . 

E 2 = — v T in (t)R' n (r)e m *—P e (cos8) B 2 = — _ R n (r) T ln {t) Pi (cos 9)e tm + (78) 

72/ ' ByY do y z Y 

7 777 CSC 1 17 I • ft 

Es = -TJTT?T en (t)K(r)Pi(cos9)e^ B 3 = - R n ( r )T in (t)e^ - Pi(co S 9) 

jy ByY y z Y do 



These are the electric multi-poles (except for 
£ = 0), both static and dynamic. The magnetic 
multi-poles are obtained Jl5[ similarly from ( |62| ) 
and are related to ([78]) by inserting an indepen- 



dent solution $ to ( p2| ) for Pm and performing the 
duality operation Ei — > Bi and B t — > — -Bj. 

For the sake of further works, we explicitly pro- 
vide here the magnetic multi-poles. These are: 



Bi 
B 2 
B 3 



e -^^T en (t)R n (r)Pe(cos6y m <t' 

yzy z 



im esc o 
IV'ByY 



T tn {t)R' n {r)P l { C osOy m 't' 



Ei = 

E 2 = - R n {r) T ln {t) P e (co S e)e im * (79) 

E 3 = 4v R "( r ) f in(t)e ml *-^ P e (cos9) 
y z Y do 



Of particular interest are the "fall-off properties" of 
the fields. It is noticed that the fields components 
are proportional to either R/y 2 or R'/yy'. 



A. Behavior of electromagnetic fields at early 
times 

At early times, we may consider the initial 
hyper-surface, at time t — ti, to contain sparsely 
distributed, pressure-less matter collapsing with- 
out any radiation. The emission of radiation will 
be switched on at some suitable time when the 
particles of matter collide and produce some radi- 
ation. 

Then, as initial condition, we have B = 
constant = Bq. Therefore, at t = U, we have 

,2 C ^- 1 2 B 2 e(e + i)/Y 2 . 

It is then to be noticed that the source-free elec- 
tromagnetic helds in this spacetime are to be given 
as "initial" conditions. 



B. Behavior of electromagnetic fields near 
the singular hyper-surface 

From © and @, it is clear that n c (£, t) is an 
increasing function of t in a collapse situation, ie, 
for B < and B < 0. Then, any initial mode 
with n c < n, ie, an initially exponential mode, 
becomes an oscillatory mode, n c > n, with the 
progress of gravitational collapse of matter. The 



frequency of oscillations of the fields then continues 
to increase with the progress of the collapse. In the 
limit of the singularity, ie, for t — * t s - the singular 
hyper-surface, n c — > oo and the frequency of field 
oscillations becomes infinite. 

Thus, all the I - modes (electric and magnetic, 
both), with I > 0, become oscillatory with the 
progress of the collapse irrespective of their nature 
at initial time. In particular, this is the case near 
the singular hyper-surface t — t s . Note that this 
is happening only asymptotically with time. The 
singular hyper-surface is to the infinite future of a 
co-moving observer who is also the "cosmological 
observer" for the spacetime of (||). 

This is clearly consistent with the result that 
the co-moving observer too becomes relativistic 
in the limit t — > t s . Therefore, the co-moving 
observer only sees radiation in this limit. That 
the co-moving observer becomes relativistic in the 
limit of the singular hyper-surface is irrespective of 
whether there is any black hole - the infinite red- 
shift surface - in the spacetime of @). This is also 
irrespective of whether there are electromagnetic 
fields in the spacetime of (g) or not. 

However, this is not the astrophysically interest- 
ing or relevant situation since it is reached for an 
infinite co-moving time when matter in the entire 
spacetime has attained relativistic speeds as a re- 
sult of the continued pile up of matter to the center 
of the spacetime. 
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C. Behavior of electromagnetic fields in the 
presence of a black hole 

Perhaps, what is astrophysically relevant is the 
situation of a black hole existing in the spacetime 
of (||). However, it must be emphasized that a 
black hole is only an infinite red-shift surface here 
and that it must exist in the spacetime of (Q) as 
an "initial" condition for non-singular and non- 
degenerate matter data. It is therefore likely that 
the black hole of (|J) may only be of academic in- 
terest. That is to say, the black hole of ^ may 
not be astrophysically relevant. 

Then, of some, nonetheless questionable, astro- 
physical interest is the behavior of electromagnetic 
fields when a black hole exists in (||) as an infinite 
red-shift surface. 
As seen in 



II H 2 



the necessary condition for 
the presence of a black hole in (g) is y(r) < 1. 
Therefore, we need to consider the field solutions 
(fr8|) in the range y : (y c , oo) when y c < 1. We will 
then obtain the fields within the black hole region 
for y : (y c , 1) and outside or exterior to the black 
hole region for y : (l,oo). 

In this case, the radial behavior is obtained for 
values y c < y < 1 and there is no difference in the 
temporal behavior of the fields. 



V. DISCUSSION 

In General Relativity, a continuum of "curved" 
4-dimensional spacetime geometry describes the 
evolution of matter. In its (3+l)-formulation, 
we can consider some distribution of gravitating 
matter on an "initial" spacelike hyper-surface, the 
Cauchy surface, and, from the Einstein equations, 
can obtain the temporal evolution of matter from 
that "initial" datum. 

We may select a variety of "matter datum" on 
the initial hyper-surface. For example, we may 
consider a single particle, ie, a single mass-point, 
or a blob of matter surrounded by vacuum or a 
blob of matter that is radiating, and evolve that 
matter datum using the field equations. 

Now, any such "initial" data is replaceable with 
"that" mass-point or "that" blob of matter being 
surrounded by "more" matter, this replacement 
being ad infinitum till the entire initial hyper- 
surface has matter everywhere. Note also that 
this "replacement" can be achieved in uncountably 
many different ways, ie, by distributing "more" 
matter in uncountably many ways on the initial 
hyper-surface. 

Further, each stage of the evolution is obtain- 
able from the previous stage of matter. This is the 



principle of causal development. 

Now, when the matter datum is specified over 
all of the initial hyper-surface, we obtain a "cos- 
mological" situation or spacetime. 

Different initial data could evolve to distinct 
four-dimensional spacetime geometries. Therefore, 
these spacetime geometries; of a mass-point, of a 
matter blob and the "cosmological" spacetime ob- 
tained for an ad infinitum replacement of any of 
the considered situations; possess different geomet- 
ric and, hence, physical features. 

The question then arises: Which geometric or 
physical features of these spacetimes are relevant 
to "real" objects of the observed Universe? To be 
able to make contact with physical objects embed- 
ded in the Universe, a spacetime of the object in 
question is then needed to be "cosmological" . 

Thus, a point of view can be advocated that the 
geometric features of only cosmological spacetimes 
are the ones that are relevant to actual physical 
objects of the real Universe. 

The point here is that, in General Relativity, the 
idealization of an "isolated" object comes with its 
own pitfalls of the above nature. The issue here 
is that we could, without changing the Newtonian 
law of gravitation, add two masses to produce a 
new mass in the Newtonian theory but that is not 
permissible with "spacetimes" of arbitrary nature 
in General Relativity. 

Note that the local spacetime geometry is, of 
course, Minkowskian since the cosmological space- 
time is locally flat. Moreover, it may also be 
that some features of a cosmological spacetime are 
present with non-cosmological spacetimes. 

In the above general spirit, we discussed, in this 
paper, "physical characteristics" of a radially ho- 
mothetic spacetime of (0). 

We first showed that the requirement of radial 
homothety, ie, of the existence of a radial homoth- 
etic Killing vector for a spherical spacetime, fixes 
the spacetime metric uniquely to (||). The exis- 
tence of a radial HKV allows, in accordance with 
Lie's theory, an arbitrary function of the co-moving 
radial coordinate in the metric. 

In a radially homothetic spacetime of (j^), there 
is, therefore, the maximal "physical" freedom, but 
not the geometrical freedom, compatible with the 
assumption of spherical symmetry. Thus, we con- 
centrated on the non-singular and non-degenerate 
matter data for (g). 

We then discussed the issue of the regularity 
of the center of a radially homothetic spacetime 
and discussed the behavior of physical quantities. 
The lack of regularity of the center of (|) for non- 
singular and non-degenerate data is compatible 
with the "physical" expectation that the "cosmo- 
logical" observer does not witness the formation of 
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the black hole. 

Next, we showed that a black hole of (||) is only 
an infinite red-shift surface and is a part of the 
"initial" data of the spacetime. 

Further, the curvature singularity of such a 
spacetime is a singular spacelike hyper-surface at 
which matter in the spacetime attains the speed 
of light. This spacetime singularity of (||) occurs 
to the future of every observer in (g). Hence, the 
spacetime singularity of (Q) is not visible to any 
observer. That is to say, the singularity of (||) is 
not naked locally or globally. 

These results are in agreement with the strong 
Cosmic Censorship that demands that the space- 
time singularity be not visible to any observer un- 
less and until it is actually encountered [ p~7| . 

We then discussed the steps of collapse of mat- 
ter from the initial dusty state. It was argued 
that "all" the necessary ingredients of the expected 
physical evolution of matter are obtainable in the 
radially homothetic spacetime. 

Following the Cohen-Kegeles generalization |H| 
of the flat-space Hertz-Debye potential formalism, 
we have also obtained the source-free electromag- 
netic fields, fl7q), in the spacetime of (||). In the 
limit of the singular hyper-surface, all the field 
modes are oscillatory. This behavior of the elec- 
tromagnetic fields is asymptotically reached in the 
limit of infinite co-moving time. 

In this spacetime, when self-gravity dominates 
with finality, the unstoppable collapse begins and 
leads to an eventual spacetime singularity, the sin- 
gular hyper-surface, only asymptotically. How- 
ever, such a situation is not astrophysically rele- 
vant since "all" the matter in the spacetime be- 
comes relativistic in this situation. This is ex- 
pected only in the asymptotic future. 

The radially homothetic spacetime of (^) there- 
fore provides us the spacetime of astrophysically in- 
teresting gravitational collapse problem. In many 
such collapse situations, matter may trap and 
carry radiation with it as it collapses. 

From the astrophysical point of view, the recent 
observations appear to point to the existence of a 
"null hyper-surface" in candidate objects. In par- 
ticular, [18|, |ll| have recently pointed out that po- 
tential black hole candidates and known neutron 
stars separate in two categories in log (L max ) ver- 
sus log (L m i n / L max ) plots for X-ray luminosity of 
these sources. 

The explanation for this separation is provided 
on the basis of the Advection-Dominated Accre- 
tion Flows (ADAF) that have more core luminos- 
ity in the X-rays in the case of Neutron Stars. The 
accreting matter exhibiting energy-advection en- 
counters the physical surface of the neutron star 
at which it deposits the stored energy to become 



X-ray bright. On the other hand, for a black hole, 
matter is expected to encounter no such physical 
surface and, hence, matter is not expected to be- 
come X-ray bright. 

The point is that in these models of ADAF the 
above is implemented in the form of a boundary 
condition that essentially implies that the matter 
exhibiting energy-advection accretes without en- 
countering a physical surface like that of a neutron 
star or not. 

The black hole of (||) - an infinite red-shift 
surface - is, in the co-moving frame, approached 
asymptotically by collapsing matter that crosses it 
in its rest frame. Consequently, matter falling onto 
the black hole will appear, to a co-moving observer, 
to be collapsing without encountering a physical 
surface like that of a neutron star. Hence, mat- 
ter displaying energy-advection will not deposit the 
stored energy at any surface. Such an object will, 
therefore, appear less X-ray bright as compared to 
an object with a physical surface. This is as per 
the interpretation of observations in [|l8|, [l9| . 

However, the black hole of (||) is a part of the 
"initial" data for it. Therefore, unless a black hole 
is assumed to exist, there is no reason for the non- 
existence of any hard surface of the collapsing ob- 
ject in the spacetime of (0). 

Thus, it seems that the reason(s) for the obser- 
vations used in [jl8| |l9| are, in all probability, dif- 
ferent than have been explored therein. For (|J), 
such reasons can only be explored on the basis of 
the detailed "physical" considerations leading to 
solutions to ( |34| ) and (|35|). 

In (||), any spherical "object" will collapse 
asymptotically to the singular hyper-surface. It is 
an Eternally Collapsing Object. However, collapse 
to singular hyper-surface of ^) requires the entire 
spacetime to be "relativistic" and that takes infi- 
nite co- moving time. Therefore, in (^), the "Eter- 
nally Collapsing Object" is also not relevant to the 
present astrophysical observations that have been 
used in g§ 

Then, from the results of the present work, it 
seems most likely that the astrophysically relevant 
possibility is that of "sufficiently Collapsed Ob- 
jects" accreting matter in astrophysical environ- 
ments. We are, of course, considering only spheri- 
cally symmetric situation in this paper. 

Now, from the above, it is then puzzling that no 
null membrane black holes emerged in the present 
study. The following issue, therefore, arises. 
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Issue of null membranes in the physically 
realizable collapse 

The standard scenario of gravitational collapse 
leading to a null membrane black hole is, so far, 
only an "expectation" based on plausibility argu- 
ments. The exact general relativistic spacetime of 
the gravitational collapse of a "physical" star sup- 
porting the standard scenario with a null mem- 
brane is not known. 

We have shown here that the radially homoth- 
etic, spherical spacetime of (||) describes the phys- 
ically realizable spherical collapse. 

Therefore, the standard picture of spherical col- 
lapse of matter leading to a null membrane black 
hole can be verified using (^). 

We have then shown that a null membrane does 
not arise for physically realizable gravitational col- 
lapse from any non-singular and non-degenerate 
data in (^|). However, as has been shown here, an 
infinite red-shift surface can "exist" in (0) as a part 
of non-singular and non-degenerate data. 

The standard picture of the gravitational col- 
lapse leading to a null membrane black hole can- 



not therefore be realized in this spacetime for non- 
singular and non-degenerate data. 

Then, the question is of the existence of another 
cosmological spacetime describing physically real- 
izable spherical collapse of matter and supporting 
the standard expectation of the existence of a null 
membrane black hole. 

The question will also be of the difference be- 
tween such a spacetime and a radially homothctic 
spacetime considered here. Does general relativity 
allow two inequivalent such spherical spacetimes? 
This is a fundamental issue. 
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APPENDIX A: COMPUTATIONAL DETAILS 



1. Electromagnetic fields 

Some useful Hodge dual operations on basis forms are: 



■kLU Alu 1 Alu 2 Alu 3 = 1 
*w° A lu 2 = - LU 1 A OJ 3 
kLU 1 Alu 3 = lu° A lu 2 

-k LU° A LU 2 A LU 3 = LU 1 
-k LU° — LU 1 A LU 2 A LU 3 



★ 1 = -uj° Alu 1 Alu 2 A lu 3 

-kLU A LU 3 = LU 1 A LU 2 
-kLU 2 A LU 3 = — LU° A LU 1 

* LU° A LU 3 A LU 1 = LU 2 

-kLU 1 = lu° Alu 2 A lu 3 

-k LU 3 = LU° A LU 1 A LU 2 



kLU Alu 1 = lu 2 A lu 3 

-kLU 1 A LU 2 = — LU° A LU 3 

■k LU 1 A LU 2 A LU 3 = LU° 

-k LU° A LU 1 A LU 2 = LU 3 

* LU 2 = LU° A LU 3 A LU 1 



Now, in the following we shall set Pm = in (|6^) , that is to say, we shall evaluate the electric multi-poles 
for the metric (|^). We shall also denote the potential as '5. Then, from (68), we obtain 



P = fw'Aw 1 
kP = 1< lu 2 A lu 3 = * y 2 Y 2 sin 9 dQ A d<\> 

dk P = ( * , f + 2*^ ) y 2 Y 2 sin 9 dt A dO A d(j) + ( *, r + 2*— I y l Y 2 sin 9 dr A d9 A 



(Al) 



Y 



<5 t ( 2^Y 

y 



— I lu" A lu 2 A lu 3 + 

y y 



2f 



IV'B jyBj 



lu 1 Alu 2 A lu 3 
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*d ★ P = 



'** 2*y\ ! / * r 
y y Y J \iy'B 



lyBJ 



d-kd-kP = 



y y y i j \iy o ^ya j 
v y y^ / V7y'S 7 sy 

( . V V V 2 



Y Y Y 2 y\ W 
+ 5* t + 2*, t B- + 2*5- - 2*5— + 2*5- ^- 

Y Y Y z Y I y 

W5 ryB / 



V 2/ 2/ / 7-BJ 7# 

^ + 2*A dt A # _ ( + 2 * ^ 7^, r A d9 
\ iy B jB j y Y J y 



I Y\ -yy'B 

- I * 40 + 2* - I -^dr A d</> 



d-kd* P = 



Y Y Y 2 ■ Y\ W 

5* « + 5* t + 2* t 5- + 2*5- - 2*5— + 2*5- ^ 

y y y 2 Y I y 



2/ 



* 



2*a\ csc6> 



7 5 



1 



IV'By 



-co Alo 1 



* 



2* 



y' y / 7-By^ 



i 



-w° A uj 2 



v 

* t0 + 2*,^ 



y J iByY 
Y \ esc 



uj° Auj 3 - i * t0 + 2* e 



Y \ 1 



y / y 2 Y 



a; 1 Auj 2 



Y J y 2 Y 
Similarly, 

P = *w° Auj 1 = Vyy'-fBdt A dr 



uj 1 A uj 3 



(A2) 



dP = * e 2/1/75 dO A dt A dr + * yy'^B d<j> A dt A dr = ^ uj 2 A uj° A uj 1 + CSC ^' uj 3 Auj A uj 1 



* e 3 csc6>*, 2 . 

-kdP = —— uj — — uj — * sin 

yy 



cscfl*^ d6> 



d ★ dP = sin [* # t dtAd(f> + * er dr A d(f> + (* ee + * e cot 0) d0 A d<j>] 

- esc (* : t dt AdO + * ^ r dr Ad9 + * ^ d</> A d0) 



* 



Bt , ,0 



y 2 ^ 



w u Auj + 



*, 



ly'ByY 



uj 1 A uj 3 + (* 00 + * e cot + esc 2 6**, 



?/ 2 y 2 



w z A w 



csc0* m 2 csc0*.A r x 2 
■ w Aw — ur A ur 



y 2 Y ly'ByY 



*d*dP = - [* ; 0o + cot 0*,e + esc 6 *,^] 



w° Aw 1 + 



* 



, ,0 



7 y'5yY 



w A u; 



CSC 0* j( /, r 

IV'ByY 



4>r 



w A 



csc6»*^ t , 3 *, et 1 



y2 y 



■w Aw 



y 2 ^ 



A uj 



Therefore, 
A5 = 



Y Y Y 2 .Y\ 1 

5*, u + 5* 4 + 2* it 5- + 2*5- - 2*5— + 2*5- 



y*,r 



+ 3* 



Y y 

yy"*.A 1 



Y / y 2 5 



y/2 J r^ljp.yly 



(*,ee + cot 6» *, e + esc 6> ^ M ) 



y 2 Y 2 



u)° Auj 1 
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2*, 



jBy 2 Y 



co° A to 2 



2csc6>* 3 2F* 
1 By 2 Y W AUJ ~yW 2 ~ 



g i 2 2 CSC 6»F* , 3 

w Aw t^ttt^ — w 1 A w 3 



(A3) 



Now, with the gauge term as 



G 



2* 



2F* 



we get 



jBy yY 
( 2Y^B 2Yi/B 2^Y 2 B 2^YB\ 1 



(A4) 



2*. 



2* 



yY yY 2 yY J By "f 2 B 2 y'y 

2F* 



~fBy 2 Y 



e n.2 2csc6<*.0 3 



jBy 2 Y 



y 2 Y 2 



LJ 1 A U? - 



lu^Alo 1 



2 esc 6>y# 
y 2 Y 2 



•w 1 A w 3 



(A5) 



Therefore, substituting ( |A3|) and ( |A5| ) in (|62|), it is seen that the w° Aw 1 term yields the required wave 
equation while all other terms lead to identities. The wave equation is: 



t/t/"¥, r \ y Y 



cot + CSC6** 



" / V y' ' ' r y' 2 / y' i 2B2 

and is amenable to solution by separation of variables. Dividing ( |A6| ) by '5 and setting 

* (t , r, M) = T ln (t) R n (r) P e (cos 0) e lm * 
where Pi(cos9) is an associated Legendre polynomial, we obtain 

1(1 + 1) 



Tin 


BT t 




Tin 


1 

BTt 


n 


R n 




2/" 


Rn 


\y 


y' 



Y 2 
R' 



n 
~f 2 B 2 

l\ 2 

y_ 

v 



(A6) 
(A7) 

(A8) 
(A9) 



where m, t and n are separation constants. In the above, an overhead dot denotes a time-derivative and 

an overhead prime denotes a derivative with respect to r. 
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